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ABSTRACT 

We have embarked on a survey for pulsars and fast transients using the 13-beam 
Multibeam receiver on the Parkes radio telescope. Installation of a digital backend 
' allows us to record 400 MHz of bandwidth for each beam, split into 1024 channels 

and sampled every 64 /is. Limits of the receiver package restrict us to a 340 MHz 
observing band centred at 1352 MHz. The factor of eight improvement in frequency 
resolution over previous multibeam surveys allows us to probe deeper into the Galactic 
plane for short duration signals such as the pulses from millisecond pulsars. We plan 
to survey the entire southern sky in 42641 pointings, split into low, mid and high 
^^O ' Galactic latitude regions, with integration times of 4200, 540 and 270 s respectively. 

, Simulations suggest that we will discover 400 pulsars, of which 75 will be millisecond 

pulsars. With ~ 30% of the mid-latitude survey complete, we have re-detected 223 
previously known pulsars and discovered 27 pulsars, 5 of which are millisecond pulsars. 
The newly discovered millisecond pulsars tend to have larger dispersion measures 
than those discovered in previous surveys, as expected from the improved time and 
frequency resolution of our instrument. 
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1 INTRODUCTION 

Pulsars are fascinating celestial objects, with spin periods 
ranging from ~ 1.6 ms to ~ 8.5 s and magnetic fields that 
span more than five orders of magnitude. Their study yields 
a better understanding of a variety of physics problems, from 
acceleration of particles in the ultra-strong magnetic field, 
to tests of gravity in the strong field regime. Pulsars provide 
useful probes of their environments whether inside pulsar 
wind nebulae, the centres of globular clusters or the Galac- 
tic Centre. Discovery of these unique objects has continued 
unabated for more than 40 years, and the diversity in their 
properties means that, even though nearly 2000 pulsars, 
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including ~ 170 millisecond pulsars (MSPs), are currently 
known, further discoveries are warranted. 

In the past decade alone, surveys have uncovered a num- 
ber of intriguing objects. These include the pulsar with the 
fastest rotation pe riod (in the globular clus ter Terzan 5; 
[ Hessels et al.ll2006jh the double pulsar system l|Burgav et al.l 
120031 ; iLvne et alj 12003) and it s impressive tests of general 
relativity (|Kramer et all [2006 V the 'missi ng link' between 
the lo w-mass X-ray binaries and t he MSPs (lArchibald et ail 

2009) , rotating radio transients (McLaughli n et alj l2006l h 
and young, energ etic pulsars such as PSR J1028-5819 
l|Keith et al.1 12008T ) with their subseque nt tie to gamma- 
ray p ulsars and their emission properties (|Weltevrede et all 

2010) . Furthermore, timing of MSPs p rovides highly ac- 
curate parameters (|Verbiest et al.l I2009T ) possibly leading 
to the detection of gravitational waves in the near future 
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l|Hobbs et alj|200^ ) . Finally, the discovery of a burst of radio 
emission o f unknown origin bu t seemingly at cosmologically 
distances (jLorimer et al] 120071 ) has motivated searches for 
other such events. 

Since the discovery of pulsars (jHewish et al] 1 1968T ) . 
searches have been split into targeted and untargeted sur- 
veys. Targets of interest include, for example, globular clus- 
ters, supernova remnants, the Galactic Centre, and other 
regions known to be rich in pulsars. Although targeted 
searches are highly efficient in their use of telescope time, 
it is impossible to predict the location of the majority of 
pulsars. Therefore, untargeted searches over large areas of 
sky are the only way to significantly increase the number of 
known pulsars. 

Since the early 1990s, large-area pulsar surveys with 
the Parkes telescope, such as th e low frequency Parkes 
Southern Pulsar Survey (PSPS; iManchester et al] Il996h 
and the hig her frequency Parkes M ultibeam Pulsar Sur- 
vey (PMPS; IManchester et a.1.1 l200lh. Swinburne Interme- 
diate Latit ude Survey f SILS;lEdwards et aljlioOlh and its 
extensions (|jacobv et a.1.1 120091 ') have led to the discovery of 
nearly 1000 pulsa rs. The all-sky PS PS was highly successful 
m finding MSPs ^ ne et alj [l998), but the low observing 
frequency limited the survey volume, especially towards the 
Galactic plane. The success of the later surveys was due to 
a combination of the higher observing frequency, the low 
system temperature and the large fi eld of view provided 
by th e 13-beam 'Multibeam' receiver ()Stavelev-Smith et al] 
Il996h . These surveys had relatively coarse frequency reso- 
lution however, and this limited the searchable volume for 
short duration pulses from MSPs and/or transient bursts. 
For the PMPS this limit was due to the costs associated with 
replicating complex sets of analogue filters over 13 beams 
and the processing power available in the late 1990s. 

In the 14 years since the commissioning of the Multi- 
beam receiver, there has been a great advance in the de- 
velopment of backend recording technology. As pen chart 
recorders gave way to digital data analysis in the early 1970s, 
so analogue signal processing equipment is now being re- 
placed by digital systems that offer a more affordable way to 
obtain higher time and frequency resolution. Increased fre- 
quency resolution allows for better correction of the disper- 
sive effects of the interstellar medium, which, in combination 
with the higher time resolution, gives increased sensitivity 
to short duration pulses, such as those from MSPs. 

Encouraged by these developments, we conceived the 
'High Time Resolution Universe' (HTRU) survey, motivated 
by two main drivers: to increase our understanding of the 
population of MSPs and to characterise the transient sky 
on timescales down to tens of microseconds. These drivers 
led us to build a new digital backend system connected to 
the Parkes Multibeam receiver. We use this new backend 
to give the HTRU survey four times the time resolution 
and eight times the frequency resolution of the PMPS. This 
allows us to detect MSPs for which interstellar dispersion 
had hidden from previous surveys (Bates et al. in prep). 
With the discover y of the radio magnetar PSR J1622-4950 
l|Levin et al] 120101 ), we have also shown that the transient 
nature of some exotic pulsars gives us good reason to re- 
survey areas of sky that are already well studied. 

The HTRU survey intends to be an all sky survey for 
pulsars and short duration radio transients, with a strong 



Table 1. Survey parameters for the three parts of the HTRU 
survey. 



Survey 


Hitrh 


Mid 






<5 < +10° 


-120° < I 


-80° < I 


Region 




I < 30° 


I < 30° 






|6| < 15° 


|6| < 3.5° 


Tobs (s) 


270 


540 


4300 


-^b earns 


443287 


95056 


15990 


7"samp (A^s) 


64 


64 


64 


B (MHz) 


340 


340 


340 


A^chan (kHz) 


390.625 


390.625 


390.625 




870 


870 


870 


Data length (samples) 


~2 22 


~2 23 


~2 26 ' 


Data/beam (GB) 


1.0 


2.0 


16.0 


Data/total (TB) 


435 


190 


250 



focus on the lower Galactic latitudes, where we make most 
use of the higher frequency resolution for negation of in- 
terstellar dispersion. A parallel effort is taking place at the 
Effelsberg radio telescope, covering the northern sky to a 
similar sensitivity, so as to provide a true all-sky survey at 
high time resolution. 

In Section 2 of this paper we describe the survey strat- 
egy and the predicted results. Section 3 outlines the observ- 
ing setup in hardware and software. Section 4 describes the 
analysis software that is being used to process the survey 
data. In Section 5 we compute and measure the effective 
sensitivity limits of the survey. Finally, in Section 6 we give 
the basic parameters of the first pulsars discovered in the 
HTRU survey. 



2 SIMULATION AND SURVEY STRATEGY 

The HTRU survey of the southern sky is composed of three 
parts, outlined in Table [1] Other than the observation time 
and sky coverage, the observing parameters are identical for 
all three parts. The low-latitude component covers a thin 
strip of the Galactic plane, with Galactic latitude |6| < 3.5° 
and longitude -80° < I < 30°. The 4200 s observation time 
is twice that of the PMPS, providing the most thorough sur- 
vey of the inner Galactic plane to date. The mid-latitude 
component covers |6| < 15° and —120° < I < 30°. Here the 
observation time is 540 s, allowing for a large area of sky to 
be covered quickly. Finally, the high-latitude component 
covers the entire sky south of declination +10° with 270 s in- 
tegrations, excluding the region covered by the mid-latitude 
survey. Broadly speaking, the aims of the low-latitude part is 
to discover faint pulsars deep in the Galactic plane, whereas 
the mid-latitude survey will find bright MSPs suitable for 
timing array projects. Finally, the high-latitude survey will 
give us a snapshot of the transient sky at 64 fis resolution. 

Detection of MSPs is limited by a combination of dis- 
persion measure (DM) broadening in individual frequency 
channels, scattering in the interstellar medium and lumi- 
nosity. For sufficiently luminous MSPs, the PMPS was pre- 
dominantly limited by dispersion broadening in the 3 MHz 
channels, rather than by scattering. The increased frequency 
resolution of the HTRU survey means that the broadening 
is now dominated by scattering in essentially all directions. 
Figure[T]shows contours of constant pulse broadening for the 
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Table 2. Simulated results for the mean number of total pulsar 
detections and new discoveries for three previous surveys and the 
three components of the HTRU survey. The numbers are split 
the simulation of 'normal' pulsars and 'MSPs' as described in 
the text. The numbers in parentheses are the actual numbers of 
pulsars found in each of the complete surveys, where we defined 
a pulsar as an MSP if it had a period of less than 30 ms and a 
period derivative less than 10 — 17 . 



Survey 



Normal 



MSPs 





Total 


New 


Total 


New 


PSPS 


281 (279) 


-(89) 


20 (19) 


-(17) 


PMPS 


1008 (1047) 


876 (721) 


19 (20) 


18 (21) 


SILS 


137 (158) 


57 (61) 


10 (11) 


6(8) 


HTRU 










Low-lat 


957 


260 


51 


33 


Med-lat 


831 


52 


48 


28 


High-lat 


783 


11 


65 


13 



Distance (kpc) 



Figure 1. Contours of constant pulse broadening timescale for 
lines of sight at b = for a frequency of 1352 MHz. The val- 
ues for the PMPS are shown in dashed lines and for the HTRU 
survey in solid lines, with contours at 0.5, 1, 5 and 10 ms. Scat- 
tering timescales and DMs were computed using the model of 
ICordes fc Laziol j2002T) . The Earth and the approximate location 
of the Galactic Centre (GC) are shown with crosses. 

PMPS and HTRU surveys, taking into account both DM 
and scattering. Even in the worst case, towards the Galactic 
Centre, the limiting distance of the HTRU survey is greater 
than twice that of the PMPS, and significantly more in all 
other directions. This implies that the mid-latitude survey, 
with an integration time only one quarter that of the PMPS, 
should still discover pulsars that lie in previously surveyed 
regions. The deeper low-latitude survey will be able to pen- 
etrate further into the Galactic plane for both MSPs and 
longer period pulsars. 

In order to determine the discovery potential of the 
HTRU survey, we ran simulatio ns of the pulsar popu lation 
based on the model presented in lLorimer et alJ (2006) using 
the PSRPOfO software. 

For the normal pulsar population, the input parameters 
to the model were as follows: 

• A log-normal spin period distribution, with mean 
10 271 ms and sigma 10 34 ms. 

• A simple power-law for the luminosity distribution, 
with index —0.59 and low luminosity cut off at 0.1 mjy kpc 2 . 

• An exponential function for the height above the Galac- 
tic plane, with scale 0.33 kpc. 

• The radial distribution of lYusifov fc Kuciikl (|2004l ). 

• A total of 28000 pulsars, such that the simulated de- 
tection of the PMPS matched closely to the true discovery 
rate. 

The number of known Galactic MSPs is small and many 
previous surveys have been restricted by interstellar disper- 
sion, making simulation of the Galactic population of MSPs 
more u ncertain and indeed t his population was not consid- 
ered bv lLorimer et all {2006). To produce simulations con- 

1 http://psrpop.sourceforge.net 



sistent with observed detection rates, we made the following 
modifications to the model: 

• The period distribution was based on that of the known 
Galactic plane MSPs. 

• A larger scale height of 0.5 kpc was used, to bet- 
ter m atch the known MSP population {Cordes fc Chernofil 
Il997h . 

Along with the three parts of the HTRU survey, we 
also simulated the PSPS, PMPS and SILS surveys. The re- 
sults are summarised in Table [2] In total, the HTRU surveys 
should discover nearly 400 pulsars, of which 75 will be MSPs, 
which would represent a doubling of the known population 
of MSPs. Final results from the survey will also allow us to 
obtain a better description of the radial and velocity distri- 
butions of MSPs. 



3 OBSERVING SYSTEM 

Observational data are acquired using the Pa rkes 21-cm 
Multibeam receiver (|Stavelev-Smith et alj[l99rj ) and a new 
digital backend system, the details of which are described in 
the following sections. 

3.1 Analogue systems 

The Parkes Multibeam receiver consists of 13 feeds at the 
prime focus of the Parkes 64 m antenna, organised as a cen- 
tral feed surrounded by two hexagonal rings, with beam 
pattern on the sky as shown in Figure [2] Our measure- 
ment of the system temperature of the central beam - 23 K 
at 1400 MHz - is slight ly higher than the value quoted in 
iManchester et all UpOl), which may be due to the refur- 
bishment of the analogue electronics in 2006 (J. Reynolds 
2010, private communication). The central feed has a sym- 
metric beam pattern, and the feeds further from the prime 
focus have a slight ellipticity and gain degradation. The de- 
tails of the receiver and telescope setup are listed in Table 
For each beam of the receiver, both polarisations are down- 
converted from 1182-1582 MHz to 0-400 MHz. 

Unfortunately, the high-frequency part of this band is 



4 M. J. Keith et al. 




°! 

Figure 2. The beam pattern of the Multibeam receiver. The 
black circles indicate the 14' FWHM of each of the 13 beams of 
the receiver. The separation Di is also equal to this 14' FWHM, 
and D2 = v^Di The filled circles show three further pointings 
interleaved to make up the the shape from which the entire survey 
is tessellated. 



Table 3. Details of the Multibeam receiver for the central 
feed, the inner ring, and outer ring of feeds, taken from 
iManchester et al.l l|200lft . 



Beam 


Centre 


Inner 


Outer 


Telescope Gain (K Jy — 1 ) 


0.735 


0.690 


0.581 


Half-power beam width (') 


14.0 


14.1 


14.5 


Beam ellipticity 


0.0 


0.03 


0.06 


Coma lobe (dB) 


none 


17 


14 



badly affected by interference from the Thuraya 3 geosta- 
tionary communications satellite, which transmits a space- 
to-earth signal in the 1525-1559 MHz band. We have in- 
stalled low-pass hardware filters after the receiver that re- 
strict the upper end of our observing band to 1522 MHz; 
therefore, the effective observing band is 340 MHz wide, cen- 
tred on 1352 MHz 



3.2 The BPSR backend System 

The Berkeley-Parkes-Swinburne Recorder (BPSR) system 
consists of 26 digital spectrometers connected to 13 server- 
class workstations that format the data and write it to disk. 
The digital spectrometers are based on the IBOrfl platform 
developed by the CASPEE0 group at t he Universi t y of C al- 
ifornia, Berkeley, described in detail in Mc Mahonl (2008). 

Each IBOB consists of two analogue-to-digital convert- 
ers and two polyphase filters implemented on an fpgaQ. 
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Each pair samples orthogonal polarisations at 800 MHz us- 
ing 8 bits per sample, then divides each signal into 1024 
spectral channels of width 390 kHz. Each spectral chan- 
nel is detected and integrated over 25 samples to yield the 
output time resolution of 64 [is. The 32-bit output of the 
spectrometer is linearly scaled before decimation to 8 bits; 
the 8-bit data are written to a sing le UDlQ packet with a 
2048 byte payload containing 1024 channels times 2 polari- 
sations. Each packet is sent to a 10 Gb/s Ethernet (lOGbE) 
interface connected directly to a Dell PowerEdge 1950 server 
via a CX4 cable. 

On the server, data are collected into 10 s blocks and 
further processed. Both polarisation streams of each channel 
are summed and then normalised by subtracting the mean 
and dividing by the standard deviation (as measured in the 
first 10 s block of each observation). The resultant data are 
then decimated to 2 bits per sample and written to disk at 
a rate of 32 Mb/s. 



3.3 Data collection 

Each server in the BPSR cluster has sufficient disk space to 
buffer ~ 3 days of observations. For the HTRU survey, we 
have chosen to archive the data to magnetic tapes before 
removing the data from the buffer, allowing for continuous 
observing. One copy is written to tape at Parkes and an- 
other is streamed via a dedicated 1 Gb/s fibre link to the 
supercomputer located at Swinburne University. Here the 
data are also written to magnetic tape and, if there is suf- 
ficient free capacity on the supercomputer, run through the 
survey data processing pipeline in close to real-time. 

In addition to the astronomical signal, the BPSR back- 
end also records a header file and a collection of auxiliary 
data files for each beam. The header file contains informa- 
tion about the observation in an XMI0 format as described 
in Section \3. 41 The auxiliary files contain bandpass and low 
resolution time series snapshots, as described in Section 13.51 
The raw time series is stored in the SlGPROcQ 'filterbank' 
data format, which consists of a binary header followed by 
continuous raw data. 

For the HTRU survey, all data are recorded using two 
bits of precision and with four samples packed per byte, with 
the first sample in the lowest significant bits of the byte. The 
data are ordered first by frequency then by time, each time 
sample being stored as a consecutive 256 byte word. The raw 
data are preceded by a 345 byte header, used for backwards 
compatibility with the SIGPROC software, however more in- 
formation is available in the XML header file. These data 
are then packed with the header and auxiliary files into a 
single GNU tarball per beam. 



3.4 Psrxml header format 

The psrxml header is a format-independent XML header 
file that describes both the observational properties of the 
recorded data and the format of the actual data itself. The 



5 User Datagram Protocol 
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psrxml header file is stored separately to the raw data, al- 
lowing trivial indexing of header data for very large data 
files, such as those produced by the BPSR backend. 

The psrxml data format is intended to be suitable for 
long term archival of pulsar data for future uses. Features 
of note are: 

• Checksumming of data blocks: As the data are writ- 
ten to disk, each block of data (usually a few hundred 
megabytes) is passed through the SHA-l|f| hash algorithm 
and the hash value stored in the header file. This allows for 
the data to be checked for errors and consistency at a later 
time. Archival data are valuable: This check-summing en- 
sures that data extracted from an archive are identical with 
the data recorded at the telescope. 

• Format and machine independence. The psrxml header 
file is written in a well-defined format, typically encoded in 
standard 7-bit ASClfl This means that the header will be 
readable on almost any computer that exists or is likely to 
exist. Since the header completely describes the data format, 
in terms of bit and byte ordering, and which bytes to read, 
one can therefore read the raw data on any computer with 
the appropriate software. 

• Human readable. XML documents can be easily read 
using a standard text editor. With the use of XSlQ 

view- 
ing a psrxml header file in a web browser displays the data 
in a pleasant format without the use of any specialist soft- 
ware. This does not change the file in any way, so the file is 
readable in ideal formats by both machine and human. 

• Extendable. Since the header is written in XML, it is 
trivial to add new XML tags to the file in a safe manner, 
storing either fields specific to a certain project, or extending 
the data format to new generic uses. 



3.5 Auxiliary files 

In addition to the header and data files, data generated for 
on-line monitoring of the system are stored for future use. 
This consists of the passband and time series measured in 
the pre-normalised, dual-polarisation, 8-bit data samples, 
completely averaged in frequency. 



3.6 Data tracking 

At the completion of every observation, a copy of the psrxml 
header file is kept on disk and indexed in a book keeping 
database. The database links the observation against the 
planned observations and marks that observation as com- 
pleted, removing that pointing from the pool of pointings 
from which observing schedules are selected. 



4 PROCESSING 

Processing of the survey is currently being carried out at 
the Jodrell Bank Centre for Astrophysics and the Swinburne 
Centre for Astrophysics and Supercomputing. 

The analysis of our data is limited by the available com- 
putational power, therefore re-processing of the data with fu- 
ture computer hardware will likely yield additional results. 
We follow a typical pulsar search work flow as described, for 
example, in iLorimer fc Kramerl |2005l ). 

4.1 Software Pipeline 

For our survey processing we have developed the hitrun 
processing pipeline, hitrun is designed to be a flexible sys- 
tem that can process both BPSR data and legacy 1-bit data 
formats without modifying the software. This flexibility, as 
well as being a useful tool, allows us to easily test the soft- 
ware on legacy data and against legacy codes. 

The hitrun pipeline is composed of several stages: Ra- 
dio frequency interference (RFI) removal, dedispersion, time 
series analysis, candidate sorting, candidate folding and op- 
timisation. Here, we outline the algorithms used for each 
step of the process. 

4-1.1 RFI removal 

RFI removal consists of two tasks, removal of 'bad' spectral 
channels, and removal of 'bad' time samples. 

First, we remove spectral channels that are affected by 
periodic RFI. Initially, we have a list of channels that are 
known to have poor data quality, and these are always re- 
moved. Then, we Fourier analyse each of the channels in- 
dividually, using the Fourier time series analysis described 
below. If any signals are detected with a signal-to-noise ratio 
(SNR) of greater than 15, we reject that channel. 

Next, we sum the data in frequency, producing a single 
time series at DM = 0. This time series is then searched 
for bursts of emission, and if found, those samples are re- 
moved from the time series. To remove longer duration RFI 
bursts, we sum together 2, 4, 8 16 and 32 adjacent samples 
and repeat the search. Any time samples that are removed 
in this process are then removed from the original data file, 
replaced by samples that are drawn from a random popu- 
lation statistically indistinguishable from the data. An ex- 
ample of this RFI removal is shown in Figure [3] where a se- 
ries of short-duration wide-band signals have been removed 
from the data. It should be noted that this process will re- 
move low-DM astrophysical signals of sufficient intensity and 
width. With our 5-cr threshold for removal, astrophysical 
signals of maximum recordable intensity will be removed if 
their DM is less than 0.12(W764^s) cm" 3 pc, where W is 
the approximate width of the pulse. 



8 Secure Hash Standard, Federal Information Processing Stan- 
dards Publications, 180-1, 1995 
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4-1.2 Dedispersion 

Since the DM of any signal is not known a priori, we have 
to search many trial DM values. We must therefore trans- 
form the single time series with many frequency channels 
into many frequency-summed time series with different dis- 
persion delays. For our survey, the processing time is largely 
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Figure 3. The intensity of recorded signals over a short time as a 
function of frequency, before (left) and after (right) RFI removal. 
The RFI can be seen as a series of vertical stripes, which are 
replaced with random noise by the RFI excision routine. 

dominated by the dedispersion time, making it important to 
choose the most efficient algorithm possible. 

The most efficient dedispersion al gorithm know n is the 
'tree' dedispersion approximation of iTavlorl (ll974T ). This, 
however, assumes a linear dependence of delay on fre- 
quency, which is a good approximation for small fractional 
bandwidths. For many recent surveys with large fractional 
bandwidths this approximation is inappropriate, and work- 
arounds have been used. In the PSPS, a two-step approach 
was used where smaller sub-bands were dedispersed us- 
ing the tree algorithm, before a final dedispersi on step to 
combine the sub-bands (Ma nchester et all 11996), In addi- 
tion to this, the PMPS used a linearisation approach where 
the data were padded with empty channels in a manner 
such that the dispersio n sweep could be considered linear 
iManchester et al.1 (|200ll ). For the HTRU survey we have de- 
cided not to use the tree algorithm as the cost of linearising 
the data outweighed the gain from the more efficient algo- 
rithm. Instead, we have chosen to implement an efficient v~ 2 
dedispersion algorithm. 

Much of the cost of dedispersion is due to the large vol- 
ume of data that must pass through the system, therefore we 
compact the data as much as possible. Since we sum 1024 2- 
bit values per output sample, there are 4096 possible output 
levels. Our processing systems perform 64-bit integer arith- 
metic, so we can pack 4 samples per 64-bit word, giving a 
factor of 4 processing throughput increase. By dividing the 
output by 16, the number of levels is reduced to 256, with lit- 
tle sensitivity loss. We therefore only need to write out 8-bits 
per sample, keeping the input/output costs to a minimum. 
Additionally, we make use of modern multi-core processors 
with multi-threaded code performing 4 dedispersion threads 
simultaneously. In future, it will likely be possible to gain 
additional performance benefits with technologies such as 
graphics processing unit co-processing. 

4-1.3 Time series analysis 

By this stage of the processing we have a number of dedis- 
persed time series on disk, which have to be searched 
for pulsar signals, hitrun provides both time-domain and 
frequency-domain methods. 

In the frequency domain we take the power spectrum 
of the data using a fast Fourier transform. We then remove 
the strong red-noise component of the spectrum at frequen- 
cies below about 10 Hz by subtraction of a running mean 
and division by a running variance. This also has the effect 



of normalising the spectrum. The resulting power spectrum 
is then searched for signals with a signal-to-noise ratio of 
greater than 6. To reclaim power in the harmonics of nar- 
row pulses, we also search spectra which have had 2, 4, 8 
and 16 harmonics summed. Potential signals from each DM 
trial are written out for later sorting, as described below. 

In the time-domain we cond uct a search with simi- 
lar me thodology to that reported in lBurke-Spolaor fc B ailcs 
(2010). In brief, we search for bright, single pulses by col- 
lating all points above 6-cr from the mean, after which the 
detections from a single pointing are compared across all the 
beams to identify sky-localised events. Details of the single- 
pulse-specific analysis and the results of the time-domain 
search will be reported in a future paper. 

4-1-4 Candidate sorting 

Candidate sorting is the process of collecting the Fourier 
detections from each DM step and combining those into a 
manageable number of candidate pulsars for further analy- 
sis. This is done by grouping detections with the same fre- 
quency into a single candidate, and by combining candidates 
that are harmonically related. In hitrun, this is done using 
software from pulsarhuntefPI 

4-1-5 Candidate folding and optimisation 

We then fold the data, producing an averaged pulse profile 
for each candidate. Folding is a coherent technique unlike 
the incoherent harmonic summing in the Fourier domain. 
HITRUN utilises the dspsr^I software to fold each candidate, 
producing profiles for each of 16 frequency bands and 32 

sub-integrations. 

Then the psrchive package l|Hotan et al.ll2004r ) is used 
to sum the folded sub-integrations and periods with a series 
of trial periods and DMs around the values returned from 
the Fourier domain search. The optimised period and DM, 
i.e. that which give highest folded SNR, are then used to 
make the final diagnostic plots for the candidate. 

4.2 Selection of candidates 

The goal of the data processing is to select a number of 
likely pulsar signals for re-observation. To achieve this, we 
take the optimised candidates and the associated diagnos- 
tic plots and view them with graphical ca ndidate selection 
tools, such as JReaper ( Keith et ai1l2009l ). These tools re- 
duce the large number of candidates to a series of points on 
a scatter chart, usually with period and SNR on the axes. 
Selecting these points presents the user with the full can- 
didate details and diagnostic plots, at which point the user 
must exercise judgement on whether the candidate should 
be flagged for further follow up. 

Since manual candidate selection is a time consum- 
ing and potentially error prone process, we are investigat- 
ing the potential for using artificial neural network (ANN) 
algorithms to pre-se lect promising candidates for viewing 
l|Eatough et al.ll2010l ). An ANN is a mathematical construct, 

11 http:/ /www. pulsarastronomy.net/wiki/Software/PulsarHuntcr 

12 http://dspsr.sourceforge.net 
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consisting of a set of nodes that are connected to form a de- 
cision tree that, based on a given input vector, derives a 
conclusion about the validity of a hypothesis. In our case, 
we provide the ANN with a vector of 'scores', formed from 
analysis of the diagnostic plots, e.g. the agreement of the ob- 
served SNR-DM curve of a candidate with theoretical pre- 
dictions as measured by a x 2 estimator. The more types of 
scores that are used, the larger (usually) the chance to dif- 
ferentiate between pulsar and non-pulsar signals, however 
the m ore known pulsar si gnals are required to 'train' the 
ANN jEatough et alj|20ich . Indeed, our ANN is poor at se- 
lecting binary or millisecond pulsars, because of the limited 
size of training set available for these classes of pulsar. The 
effectiveness of this technique can be seen in the discovery 
of PSR J1701-44, which was identified as a candidate by 
the ANN. The SNR of the candidate was rather low, at only 
7, and would quite probably have been missed by a graph- 
ical interface reliant upon a human to select the candidate 
manually. 

Although the ANN is limited in the types of pulsars 
it will select, it does provide a thorough analysis of all the 
candidates free from human error. We find that adding the 
ANN to our existing candidate selection effort gives a ben- 
eficial additional way to select candidates, with little cost 
in terms of required computing power. Direct inspection of 
candidates is vital, however, and remains the primary way 
that pulsars have been discovered in the HTRU survey. 



5 SURVEY SENSITIVITY 

Knowledge of the minimum detectable flux density of the 
survey is vital for determining the survey effectiveness. De- 
termining this limiting flux density is not trivial, since the 
survey is not uniformly sensitive, even over a small region 
of sky, nor to all pulsars. One option is to compute a theo- 
retical mean sensitivity to some standard 'typical' pulsar, as 
done by other large pulsar surveys. This gives a value that 
can be easily compared to other similar surveys. 

To measure the actual sensitivity of the survey, we can 
look at re-detections of known pulsars. Since published val- 
ues exist for the position and flux density of thousands of 
pulsars, we can cross match those positions with our survey 
pointings. We then compute the theoretical SNR value for 
that pulsar given the offset in the beam position and com- 
pare that to the output of the search pipeline. We consider 
each of these in in turn below. 



5.1 Theoretical sensitivity 

The radiometer equation can be used to compute the fun- 
damental limiting flux density for the centre of the central 
beam of a given survey observation. 



(1) 



where a is the minimum acceptable SNR, T sys is the system 
noise temperature, T s k y the sky noise temperature, G is the 
system gain, B is the observing bandwidth and T b s is the 
integration time. T s k y is only a significant contribution to 
the system noise for low Galactic latitudes and has a strong 
dependence on the area of sky that is observed. The value of 



Table 4. The maximum, minimum and mean sky noise temper- 
ature over the three survey regions. Note that the values for the 
mid-latitude survey exclude the low-latitude survey region. 



Survey 


Mean (K) 


Max (K) 


Min 


Low 


7.6 


36.0 


1.6 


Mid 


2.5 


9.1 


0.6 


High 


1.0 


2.1 


0.6 



T a k y at 1352 MHz for the survey region is extrapolated from 
the lHaslam et ail (|l98lh sky map at 408 MHz. For Equa- 
tion [T] we use the mean values listed in Table [4] however it 
should be noted that T s k y , and therefore S m i n , varies consid- 
erably across the sky. Here we use a — 8, T sys = 23 K and 
G = 0.735 to compute a S min of 0.20, 0.47 and 0.61 mjy 
for the low, mid and high latitude portions of the survey 
respectively. However as described below, the practical lim- 
iting flux density is lower due to the narrow pulse width of 
most pulsars. 

Our Fourier based detection algorithm is not uniformly 
sensitive to all DMs, pulse periods and pulse shapes. This is 
somewhat difficult to quant i fy, how ever we follow the ap- 
proach of Manches ter et ail (|200ll ) by modelling the fre- 
quency response of the pulsar as a uniform train of impulses 
separated by frequency f p = 1/P. This train of pulses has 
amplitude 1/S m m, to which we multiply a series of func- 
tions representing the various transformations that the data 
undertake. 

We can model each pulse as a Gaussian with width 
W50 = 0.05P, so that the frequency response is multiplied 
by the Gaussian 



F(f) = exp 



(7V.fW 50 ) 2 

41n(2) 



(2) 



The pulses are also convolved with the effect of interstel- 
lar dispersion, which can be modelled as another Gaussian 
with width tdm, the dispersion delay. Finally we must ac- 
count for the finite sampling interval, r s , which multiplies 
the amplitudes by 



F(f) 



sin (nfr 3 



tt/t s 



(3) 



The final SNR is computed after summing N harmonics, 
where N € {1,2,4,8,16}. Figure g] shows the minimum de- 
tectable flux density as a function of pulsar period and DM 
for the mid-latitude survey. It can be seen that the HTRU 
survey is sensitive to MSPs out to DMs of a few hundred 
cm _3 pc. The effects of T a k y and the differing observing time 
for the three surveys can be seen in Figure [5] which shows 
the sensitivity curves for the three surveys at the minimum 
and maximum T s k y in that survey region. 

The mean sensitivity of the survey is reduced by a num- 
ber of factors. Firstly, the outer beams of the receiver are 
less sensitive than the centre beam, and the sensitivity of 
each beam reduces towards the edges. This means that the 
sky coverage is not uniform, with the mean sensitivity over 
the mid-latitude survey of about 0.25 mjy. 
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Figure 4. The mean limiting flux density for the mid-latitude 
survey shown as a function of pulsar period for various values of 
DM. Values shown here assume a pulse duty cycle of 5%. 
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Figure 5. Comparison of the sensitivity of the three survey 
regions for a pulsar with a DM of 100 cm - 3 pc. Each region has 
two curves, showing the values for the minimum and maximum 
value of T s ]j y in the survey region. 

5.2 Re-detection of known pulsars 

A total of 290 previously known pulsars lie within one 
beamwidth of one of the HTRU survey pointings collected to 
date. For each pulsar we obtained the flux density and its er- 
ror from the pulsar cataloguj^l. a value of T s k y at the pulsar 
position, and the offset between its position and the centre 
of the relevant beam. We then computed the expected SNR 
for each pulsar following the process outlined in Section \5. II 
In all, we detected 223 of the 290 pulsars. For the 67 
pulsars we failed to detect, 5 observations were ruined by 
broad-band RFI, 54 pulsars were not detected above a SNR 
of 10 even after we folded observations at the known pul- 
sar ephemeris. This leaves 8 pulsars that are seen in folded 
observations with a SNR of greater than 10 but were not de- 
tected in our search procedure. Each of the the 8 pulsars had 
a spin period above 500 ms, and so the Fourier SNR may 
well have been reduced by a larger red-noise component in 
the data, or the pulsar period falling between Fourier bins. 

13 http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ 
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Figure 6. The predicted SNR for 223 known pulsars re-detected 
in the HTRU survey, shown as a function of the observed SNR. 
Errors shown are from the published error in the pulsar flux den- 
sity. The line shows the expected 1:1 correlation. 

The loss of pulsar detections due to RFI is an uncontrollable 
aspect of search observations. Although we have applied the 
techniques described in Section l4.1.1l they are unable to re- 
move this type of br oad-band RFI. Future re-processing us- 
ing the technique of lEatough et all (|2009h may be possible, 
however this is currently unfeasible as it would increase the 
time for processing the survey data by more than a factor 
of ten. 

Figure |S] shows the expected versus actual SNR for the 
223 detected pulsars. Although the correlation between the 
predicted and observed values is good, the predicted values 
are, on average, somewhat higher than the observed values. 
This can be caused by either sensitivity loss in our observ- 
ing system, or error in the catalogued flux densities. First, 
we checked that the backend was not introducing any loss 
by collecting data simultaneously with BPSR and the Parkes 
Digital Filterbank 3 (PDFB3). These observations show that 
the BPSR and PDFB3 performance are equivalent and we 
therefore consider it unlikely that the backend is introduc- 
ing any losses. Hobbs (private communication), has shown 
that the catalogue flux densities appear to be systematically 
higher than those obtained in interferometric observations. 
This bias may be due to the tendency to publish only the 
data with the highest SNR, for example when scintillation 
boosts the apparent flux density. Given our ability to de- 
tect the majority of the pulsars that are in our search area, 
we are relatively confident that the system is performing to 
specification. 



6 NEW PULSAR DISCOVERIES 

In this paper we present the results of the analysis of the 
first ~ 30% of the mid-latitude part of the survey. Results 
so far include: 
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Table 5. Basic parameters for the 27 pulsars discovered in the 
HTRU survey to date. Pulsars for which no full timing solution is 
published have been assigned a temporary name containing only 
two digits of declination. 



PSR J 


I 

(°) 




b 

(°) 




Period 
(s) 


DM 

(cm~ 3 pc) 


1 1 OR — RR 1 

1 1 Z. J JU 


291 


8 


2 


6 





003102 


125 


1330-52 


308 


8 


9 


6 





648104 


149 


1433-50 


318 


9 


9 


1 


1 


017495 


98 


1442 — 51 


320 


1 


7 


7 





732061 


97 


1517—46 


327.38 


9.22 





88661 


139 


1 fi1 9 — RR 

1 >J L Z, JO 


327 





-5 








615520 


172 


1623-4950 2 


333 


9 


-0 


1 


4 


325800 


881 


1625-49 


334 


g 


o 


o 





355856 


717 


1 697— 59 


327 


3 


-7 


3 





354239 


92 


1 fi4.fi — 3fi 


347 


1 


5 


8 





212316 


222 


1 fisn— fin 


328 


2 


-10 


2 





583771 


106 


1701-44 


342 


3 


-1 


3 





755535 


413 


1705-43 


343 


3 


-1 


5 





222561 


184 


1705—61 


328 


7 


-12 


2 





808546 


87 


1708-35 1 


350 


4 


3 


2 





004505 


146 


1731-18 1 


6 


9 


8 


2 





002345 


106 


1753-38 


352 


3 


-6 


4 





666804 


168 


1754-24 


4 


9 





8 


2 


090263 


789 


1801-32 1 


358 


9 


-4 


5 





007454 


177 


1802-33 


357 


6 


-5 


5 


2 


461052 


254 


1803-33 


358 





-5 


6 





633412 


171 


1810-01 


27 





8 


6 





744976 


135 


1811-24 1 


7 


1 


-2 


5 





002661 


60 


1811-49 


344 


2 


-14 


2 


1 


432704 


49 


1812-27 


4 





-4 


6 





236983 


105 


1814-05 


23 


8 


5 


8 


1 


014405 


119 


1854-15 


19 





-8 





3 


453121 


157 



1 Timing parameters in Bates et al. (in prep). 

2 Timing parameters in Ijevine^il^ (2010). 



• Discovery of 5 MSPs, all of which are in binary systems. 
Full details of these pulsars will be reported in Bates et al. 
(in prep). 

• Discovery of the radio pulsar with the highest known 
magnetic field which appears to be a radio-loud magnetar, 
as reported in iLevin et al.1 (|20ld ). 

• Discovery of 21 pulsars with periods between 0.2 and 
3.5 s and moderate to high DMs. These are likely 'standard' 
solitary radio pulsars but full timing solutions still need to 
be obtained. 

The Galactic coordinates, spin period and DM for each of 
these pulsars is given in Table [5] 

6.1 Comparison to PMPS 

Figure [7] shows a comparison between the sensitivities of the 
mid-latitude and low-latitude HTRU survey and the PMPS 
for a DM of 100 cm" 3 pc. Although the mid-latitude HTRU 
survey has only one quarter the observing time of the PMPS, 
it is significantly more sensitive to short period pulsars at 



\ PMPS 

0.8 - 




HTRU: Low-latitude 



1 10 100 

Period (ms) 

Figure 7. Comparison of the Mid and Low latitude parts of 
the HTRU survey and the PMPS. The lines give the sensitivity 
limit for a DM of 100 cm - 3 pc. The shaded areas show the range 
of sensitivities for DMs between 50 and 150 cm - 3 pc. Note that 
the limiting flux density for the Mid-latitude survey used for this 
curve is computed using the mean T s ^ y for the PMPS survey 
area. Also shown as diamonds, are the periods and estimated flux 
densities for the four MSPs discovered in HTRU survey that lie 
in the PMPS region. 

large DMs. This is due to a combination of the faster sam- 
pling rate, narrower frequency channels, slightly increased 
total bandwidth and multi-bit sampling of our survey. 

Of the new discoveries, 9 pulsars lie within the sky cov- 
erage of the PMPS. We have re-analysed archival data for 
each of these pulsars and the results are presented in Ta- 
ble [6] Four of these pulsars are MSPs with high DMs and 
are included on Figure [7] It can be seen that all bar one 
fall below the (nominal) PMPS detection threshold. In fact, 
re-examination of the PMPS archival data shows that both 
PSR J1708-35 and J1801-32 are visible but that the pulse 
broadening from the 3 MHz channel bandwidth makes them 
difficult to distinguish from sinusoidal interference in a blind 
search. 

Of the remaining five pulsars, three are seen in the 
archival data with low SNR and none were selected for re- 
observation in the PMPS. Although the pulsar positions are 
not well known at this time, the nominal position of each 
falls more than 5 arc-minutes from the beam centre of the 
PMPS observation. Therefore we suggest that these pulsars 
fall in the regions between beam centres where the PMPS 
sensitivity is lowest. Two pulsars appear to be below the 
detection limit of the PMPS including the newly discovered 
magnet ar PSR J1622-495 0, which shows extreme time vari- 
ability (jLevin et al.ll2010t ). 



7 CONCLUSION 

We have described the system configuration and initial dis- 
coveries of a new pulsar survey using the Parkes radio tele- 
scope. At the time of writing, a total of 250 pulsars have 
been detected, of which 27 are new discoveries. The high 
time and frequency resolution of our digital backend sys- 
tem leads to increased sensitivity to short-period, high-DM 
pulsars compared to previous surveys. Indeed, 5 of our new 
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Table 6. Pulsars discovered in the HTRU survey for which 
archival data exists in the PMPS at the position of the pulsar. The 
uniqu e PMPS grid identification number (see [Mancheste r et al.l 
l200ll for details) for the processed beam, the MJD of the obser- 
vation and the SNR and periodicity if the pulsar was detected. 



PSR J 




Grid ID 


MJD 




SNR 


Period (s) 


1125- 


58 


G4707513 


51296 


.42 










1622- 


4950 


G4887499 


50686 


.25 










1625- 


49 


G4891500 


50686. 


.23 


7 


.9 





.3558530 


1701- 


-44 


G4924493 


51152. 


,94 


6 


.4 





.7555360 


1705- 


-43 


G4928493 


51457. 


,10 


10 


.4 





.2225614 


1708- 


-35 


G4959516 


51718. 


,68 


17 


.4 





.0045052 


1754- 


-24 


G5021504 


50843 


,07 










1801- 


-32 


G4995477 


52080 


,69 


14 


.4 





.0075344 


1811- 


24 


G5030488 


51245 


,84 











discoveries are MSPs at high DM, all of which lie within 
previously surveyed regions of sky. 
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